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Abstract 

The detection of high-energy particles, cosmic rays (CRs), deep inside 
the heliosphere implies that there are, at least, three distinctly different 
stages in the lifetime of a CR particle: acceleration, propagation in the 
interstellar medium (ISM), and propagation in the heliosphere. Gamma 
rays produced by interactions of CRs with gas, radiation, and magnetic 
fields can be used to study their spectra in different locations. Still, ac- 
curate direct measurements of CR species inside the heliosphere (such as 
their spectra and abundances) are extremely important for the under- 
standing of their origin and propagation. In this paper, an emphasis is 
made on very recent advances and especially on those where GLAST and 
PAMELA observations can lead to further progress in our understanding 
of CRs. 
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1 Introduction 

Cosmic rays and 7-rays are intrinsically connected: 7-ray emission is a direct 
probe of proton and lepton spectra and intensities in distant locations of the 
Galaxy. Diffuse emission accounts for ^80% of the total 7-ray luminosity of 
the Milky Way and is a tracer of interactions of CR particles in the ISM pLj . On 
the other hand, direct measurements of CR species can be used to probe CR 
propagation, to derive propagation parameters, and to test various hypotheses 
[2J. Luckily, two missions of the present and near future are targeting these 
issues and present unique opportunities for breakthroughs. The Payload for 
Antimatter-Matter Exploration and Light-nuclei Astrophysics (PAMELA) [3] 
has been launched in June 2006 and is currently in orbit. During its projected 
3 yr lifetime it will measure light CR nuclei, antiprotons, and positrons in the 
energy range 50 MeV/n - 300 GeV/n with high precision. The Gamma-ray 
Large Area Space Telescope (GLAST) [4 is scheduled for launch in early 2008. 
It has significantly improved sensitivity, angular resolution, and much larger 
field of view than its predecessor EGRET and will provide excellent quality 
data in the energy range 20 MeV - 300 GeV. 

2 Cosmic-ray accelerators 

A supernova (SN) ~ CR connection has been discussed since the mid-1930s 
when Baade and Zwicky proposed that SNe are responsible for the observed 
CR flux. The first direct evidence of particle acceleration up to very high 
energies (VHE) came from observations of synchrotron X-rays from the super- 
nova remnant (SNR) SN1006 [5]. More recently, observations of TeV 7-rays [5] 
confirm the existence of VHE particles. Still, definitive proof that SNRs are 
accelerating protons is absent. Recent observations of the SNR RX J1713 by 
HESS suggest that its spectrum is consistent with the decay of pions produced 
in pp- interactions 'W , while the spectrum from inverse Compton scattering 
(ICS) does not seem to fit the observations. However, a calculation of ICS 
and synchrotron emission using a one-zone model [7] and a new calculation 
of the interstellar radiation field (ISRF) shows that a leptonic origin is also 
consistent with the data [8]. Another interesting case study is the composite 
SNR GO. 9+0.1 near the Galactic center [9] which can also be fitted using the 
leptonic model [8J. In this instance, the major contribution comes from ICS 
off optical photons while the 7-ray spectrum exhibits a "universal" cutoff in 
the VHE regime due to the Klein-Nishina effect. If this modelling is correct, 
GLAST observations can be used to probe the ISRF in the Galactic center. 
Observations of SNRs by GLAST will be vital in distinguishing between lep- 
tonic or hadronic scenarios as their predictions for the spectral shape in the 
GeV energy range are distinctly different. 
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A new calculation of the ISRF shows that it is more intense than previously 
thought, especially in the inner Galaxy where the optical and infrared photon 
density exceeds that of the cosmic microwave background (CMB) by a factor 
of 100 For a source in the inner Galaxy, properly accounting for inverse 
Compton energy losses flattens the electron spectrum in the source compared 
to the case of pure synchrotron energy losses [TU] ■ This effect leads to a flatter 
intrinsic 7-ray spectrum at the source. On the other hand, the intense ISRF 
also leads to 77-attenuation which starts at much lower energies than for the 
CMB alone [Tl]. For VHE 7-ray sources located in the inner Galaxy, the 
attenuation effects should be seen for energies ~30 TeV [T^ . 

A new class of VHE 7-ray sources and thus CR accelerators, close binaries, 
has been recently found by HESS [13]. The observed orbital modulation due 
to the 77-attenuation on optical photons of a companion star testifies that 
the VHE 7-ray emission is produced near the compact object. Such an effect 
has been predicted long ago in a series of papers [T3] , where the light curves 
were calculated for binaries which were suspected to be VHE emitters at that 
time, like Cyg X-3. The phase of the maximum of the emission depends on the 
eccentricity of the orbit and its orientation; for orbital parameters of LS 5039 it 
is about 0.7 [2], in agreement with observations. For a recent discussion of the 
orbital modulation in LS 5039, see [T^. GLAST observations in the GeV-TeV 
range will be the key to understanding the emission mechanism(s). 

3 Propagation of cosmic rays and diffuse gamma-ray emission 

3.1 Particle propagation near the sources 

Diffuse 7-ray emission in the TeV energy range has been recently observed 
by HESS TB] from the Galactic center. The emission clearly correlates with 
the gas column density as traced by CS. If this emission is associated with a 
relatively young SNR, say Sgr A East, observation of the individual clouds will 
tell us about CR propagation there. A simple back-of-the-envelope calculation 
shows that if the SNR age is <10 kyr and the shock speed is <10^ km/s, 
the shell size should be <100 pc, while the emission is observed from distant 
clouds up to 200 pc from the Galactic center. The emission outside the shock, 
therefore, has to be produced by protons which were accelerated by the shock 
and left it some time ago. The spectrum of such particles can be approximated 
by a (5-function in energy which depends on the SNR age; the resulting 7-ray 
spectrum is essentially flatter than expected from a power-law proton spectrum 
in the shell (Figure 1) [iTllTH]. Observations of individual clouds in the GLAST 
energy range will be a direct probe of this model and thus of proton acceleration 
in SNRs. 

Milagro has recently observed the diffuse emission at 12 TeV from the 
Cygnus region '19]. The observed emission, after subtraction of point sources 
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Figure 1: The spectrum of 7-rays [H] from the gas clouds outside of the SNR 
shell (monoenergetic protons of 25 TeV, dots) and from the shell (power-law 
with index —2.29, dashes); normalisations are arbitrary. The solid line is the 
total spectrum. Data: HESS observations of the Galactic center ridge [6]. 

correlates with gas column density. The VHE 7-ray flux is found to be larger 
than predicted by the conventional and even optimized model tuned to fit the 
GeV excess [5S] . This may imply that freshly accelerated particles interact with 
the local gas, but other possibilities such as ICS or unresolved point sources 
can not be excluded on the base of Milagro observations alone. 

These observations show that the diffuse emission exists even at VHE ener- 
gies and variations in brightness are large. A contribution of Galactic CRs to 
the diffuse emission in this energy range is still significant. Observations of the 
diffuse emission may be used to study CR propagation and their penetration 
into molecular clouds. GLAST is ideally suited to address these issues. 

3.2 Galactic cosmic rays 

Propagation in the ISM changes the initial spectra and composition of CR 
species due to spallation, energy losses, energy gain (e.g., diffusive reacceler- 
ation), and other processes (e.g., diffusion, convection). The destruction of 
primary nuclei via spallation gives rise to secondary nuclei and isotopes which 
are rare in nature (i.e., Li, Be, B), antiprotons, pions and kaons that decay 
producing secondary leptons and 7-rays. Studies of stable secondary nuclei 
(Li, Be, B, Sc, Ti, V) allow the ratio (halo size)/ (diffusion coefficient) to be 
determined and the incorporation of radioactive secondaries (^°Be, ^^Al, ^^Cl, 
^''Mn) is used to find the diffusion coefficient and the halo size separately. For 
a recent review on CR propagation see [2]- 
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Measurement of the B/C ratio with a single instrument and in a wide energy 
range is long overdue. The best data >0.8 GeV/n to-date are those taken by 
the HEAO-3 experiment more than 25 years ago, while modern spacecraft, 
e.g., ACE, provide high quality data at low energies 150-450 MeV/n [5D]. The 
sharp maximum in the B/C ratio observed at ~1 GeV/n is difficult to explain 
in a physical model and has been long debated; it may well be an instrumental 
artefact. On the other hand, the high energy tail of the B/C ratio is sensitive 
to the rigidity dependence of the diffusion coefficient and thus its accurate 
measurement can be used to distinguish between models of CR propagation in 
the ISM. The PAMELA has the capability of measuring the B/C ratio in the 
energy range 100 MeV/n - 250 GeV/n and will address both issues. 

The majority of CR antiprotons observed near the Earth are secondaries 
produced in collisions of CRs with interstellar gas. Because of the kinematics 
of this process, the spectrum of antiprotons has a unique shape distinguishing 
it from other CR species. It peaks at ^2 GeV decreasing sharply toward lower 
energies. Because of their high production threshold and the unique spectral 
shape antiprotons can be used to probe CR propagation in the ISM and the 
heliosphere, and to test the local Galactic average proton spectrum (for a dis- 
cussion and references see HJ). Because the pp (and pp) total inelastic cross 
section is ten times smaller than that of carbon, the ratio p/p can be used to 
derive the diffusion coefficient in a much larger Galactic volume than the B/C 
ratio. The CR p spectrum may also contain signatures of exotic processes, such 
as, e.g., WIMP annihilation. However, currently available data (mostly from 
BESS flights ^22J) are not accurate enough, while pubhshed estimates of the 
expected flux differ significantly; in particular, the reacceleration model under- 
produces antiprotons by a factor of ^2 at 2 GeV "51] . Secondary CR positrons 
are produced in the same interactions as antiprotons and are potentially able 
to contribute to the same topics. Accurate measurements of the CR spec- 
trum may also reveal features associated with the sources of primary positrons, 
such as pulsars. During its lifetime, PAMELA will measure CR antiprotons 
and positrons in the energy range 50 MeV - 250 GeV with high precision . 
Independent CR p measurements below '^3 GeV will be provided by the new 
BESS-Polar instrument scheduled to fly in December of 2007 [5^ . 

The diffuse emission is a tracer of interactions of CR particles in the ISM 
and is produced via ICS, bremsstrahlung, and 7r"-decay. The puzzling excess in 
the EGRET data above 1 GeV relative to that expected has shown up in all 
models that are tuned to be consistent with local nucleon and electron spectra 
|25j . The excess has shown up in all directions, not only in the Galactic plane. 
If this excess is not an instrumental artefact, it may be telling us that the CR 
intensity fluctuates in space which could be the result of the stochastic nature of 
supernova events. If this is true, the local CR spectra are not representative of 
the local Galactic average. Because of the secondary origin of CR antiprotons, 
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their intensity fluctuates less than that of protons and p measurements can be 
used instead to derive the average local intensity of CR protons. Interestingly, 
a model based on a renormalized CR proton flux (to fit antiprotons) and a 
CR electron flux (using the diffuse emission itself), the so-called optimized 
model, flts the all-sky EGRET data well providing a feasible explanation 
of the GeV excess. The GLAST observations of the diffuse emission will be 
able to resolve this puzzle. On the other hand, accurate measurements of CR 
antiprotons by PAMELA can be used to test the CR fluctuation hypothesis. 

4 Cosmic rays in the heliospiiere 

Interestingly, GLAST will be able to trace the CRs in the hcliosphere as well. 
The ICS of CR electrons off solar photons produces 7-rays with a broad dis- 
tribution on the sky contributing to a foreground that would otherwise be 
ascribed to the Galactic and extragalactic diffuse emission [25|. Observations 
by GLAST can be used to monitor the heliosphere and determine the electron 
spectrum as a function of position from distances as large as Saturn's orbit to 
close proximity of the Sun, thus enabling unique studies of solar modulation. 
A related process is the production of pion-decay 7-rays in interactions of CR 
nuclei with gas in the solar atmosphere [27] . The albedo 7-rays will be observ- 
able by GLAST providing a possibility to study CR cascade development in the 
solar atmosphere, deep atmospheric layers, and magnetic fleld(s). The original 
analysis of the EGRET data assumed that the Sun is a point source and yielded 
only an upper limit |2Sj. However, a recent re-analysis of the EGRET data [23] 
has found evidence of the albedo (pion-decay) and the extended ICS emission. 
The maximum likelihood values appear to be consistent with the predictions. 

GLAST will also be able to measure CR electrons directly. It is a very 
efficient electron detector able to operate in the range between '^20 GeV and 2 
TeV [30 . The total number of detected electrons will be ^10^ per year. Accu- 
rate measurements of the CR electron spectrum are very important for studies 
of CR propagation and diffuse 7-ray emission. There is also the possibility to 
see the features associated with the local sources of CR electrons. 

The Moon emits 7-rays [HI [SI] due to CR interactions in its rocky surface. 
Monte Carlo simulations of the albedo spectrum using the GEANT4 framework 
show that it is very steep with an effective cutoff around 3 GeV and exhibits a 
narrow pion-decay line at 67.5 MeV [32]. The albedo ffux below ~1 GeV sig- 
nificantly depends on the incident CR proton and helium spectra which change 
over the solar cycle. Therefore, it is possible to monitor the CR spectrum at 1 
AU using the albedo 7-ray flux. Simultaneous measurements of CR proton and 
helium spectra by PAMELA, and observations of the albedo 7-rays by GLAST, 
can be used to test the model predictions. Since the Moon albedo spectrum 
is well understood, it can be used as a standard candle for GLAST. Besides, 
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the predicted pion-decay line at 67.5 MeV and the steep spectrum at higher 
energies present opportunities for in orbit energy cahbration of GLAST. 

I thank Troy Porter for useful suggestions. This work was supported in part 
by a NASA APRA grant. 

References 

[1] I. V. Moskalenko, A. W. Strong, in Int. Conf. on Astrophysical Sources of 
High Energy Particles and Radiation, eds. T. Buhk et al. (NY: AIP), AIP 
Conf. Proc. 801, 57 (2005); larXiv:a stro-ph/d 509414 

I. V. Moskalenko, A. W. Strong, O. Reimer, in "Cosmic Gamma-Ray 
Sources," eds. K. S. Cheng & G. E. Romero (Dordrecht: Kluwer), ASSL 
304, 279 (2004); arXiv:astro-ph/0402243 

[2] A. W. Strong, I. V. Moskalenko, V. S. Ptuskin, Annu. Rev. Nucl. Part. 
Sci. 57, 285 (2007) 

[3] P. Picozza, et al., Astropart. Phys. 27, 296 (2006) 

[4] P. Michelson, in X-Ray and Gamma-Ray Telescopes and Instruments for 
Astronomy, eds. J. E. Triimper & H. D. Tananbaum, Proc. SPIE 4851, 
1144 (2003) 

[5] K. Koyama, et al., Nature 378, 255 (1995) 

[6] F. Aharonian, et al., A&A 449, 223 (2006) 

[7] F. Aharonian, A. M. Atoyan A&A 351, 330 (1999) 

[8] T. A. Porter, I. V. Moskalenko, A. W. Strong, ApJ 648, L29 (2006) 

[9] F. Aharonian, et al., A&A 432, L25 (2005) 

[10] J. A. Hinton, F. A. Aharonian, ApJ 657, 302 (2007) 

[11] I. V. Moskalenko, T. A. Porter, A. W. Strong, ApJ 640, L155 (2006) 

[12] T. A. Porter, I. V. Moskalenko, A. W. Strong, in 1st Int. GLAST 
S ymp., eds. S. R itz et al. (NY: AIP), AIP Conf. Proc. 921, 411 (2007); 
iarXiv:0704. 17031 

[13] F. Aharonian, et al., A&A 460, 743 (2006) 

[14] I. V. Moskalenko, S. Karakula, W. Tkaczyk, in Proc. 22nd Int. Cosmic 
Ray Conf. (Dublin), 1, 544 (1991) 

I. V. Moskalenko, S. Karakula, W. Tkaczyk, MNRAS 260, 681 (1993) 
I. V. Moskalenko, S. Karakula, ApJS 92, 481 (1994) 
I. V. Moskalenko, Space Sci. Rev. 72, 593 (1995) 



8 



Frontier Science 2007, Science with the New Generation of HE Experiments 



[15] M. Botther, C. D. Dermer, ApJ 634, LSI (2005) 
G. Dubus, A&A 451, 9 (2006) 

D. Khangulyan, F. Aharonian, V. Bosch-Ramon, MNRAS, submitted 
(2007) ; larXiv:0707. 16891 

[16] F. Aharonian, et ah. Nature 439, 695 (2006) 

[17] S. Gabici, F. A. Aharonian, ApJ Lett., submitted (2007): [arXiv:0705.30lil 

[18] I. V. Moskalenko, T. A. Porter, M. A. Malkov, P. H. Diamond, in Proc. 
30th Int. Cosmic Ray Conf. (Merida), submitted (2007): larXiv:0705.3854l 

[19] A. Abdo, et ah, ApJ 658, L33 (2007) 

[20] M. E. Wiedenbeck, et ah. Space Sci. Rev. 99, 15 (2001) 

[21] I. V. Moskalenko, A. W. Strong, J. F. Ormes, M. S. Potgieter, ApJ 565, 
280 (2002) 

I. V. Moskalenko, A. W. Strong, S. G. Mashnik, J. F. Ormes, ApJ 586, 
1050 (2003) 

[22] A. Yamamoto, et ah, Nucl. Phys. B Proc. Suppl. 166, 62 (2007) 
[23] A. M. Lioneto, A. Morselh, V. Zdravkovic, JCAP 9, #010 (2005) 
[24] S. D. Hunter, et ah, ApJ 481, 205 (1997) 

[25] A. W. Strong, I. V. Moskalenko, O. Reimer, ApJ 537, 763 (2000) 
A. W. Strong, 1. V. Moskalenko, O. Reimer, ApJ 613, 962 (2004) 

[26] L V. Moskalenko, T. A. Porter, S. W. Digel, ApJ 652, L65 (2006) 

E. Orlando, A. W. Strong, Ap&SS 309, 359 (2007) 

[27] D. Seckel, T. Stanev, T. K. Gaisser, ApJ 382, 652 (1991) 

[28] D. J. Thompson, D. L. Bertsch, D. J. Morris, R. Mukherjee, J. Geophys. 
Res. A 120, 14735 (1997) 

[29] E. Orlando, D. Retry, A. W. Strong, in 1st Int. GLAST Symp., eds. S. Ritz 
et ah (NY: AlP), AIP Conh Proc. 921, 502 (2007): [arXrv:0704.0462l 

[30] A. A. Moiseev, J. F. Ormes, 1. V. Moskalenko, in Proc. 30th Int. Cosmic 
Ray Conh (Merida), submitted (2007); arXiv:0706.0882l 
A. A. Moiseev, These Proceedings 

[31] D. J. Morris, J. Geophys. Res. A 89, 10685 (1984) 

[32] 1. V. Moskalenko, T. A. Porter, in Proc. 30th Int. Cosmic Ray Conh 
(Merida), submitted (2007): [arXi^:0705.3856l 



